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ABSTRACT 



Finite-amplitude standing waves in air at ambient con- 
ditions contained in a rigid-walled cylindrical tube with 
a large length-to-diameter ratio were experimentally inves- 
tigated. The pressure waveform at the end of the tube was 
digitized and Fourier analyzed on an IBM 360 digital 
computer. Amplitudes and phases were obtained for all 
harmonics with amplitudes greater than 1% of the fundamental 
for strength parameters from 0.25 to 1.00 and for frequency 
parameters from -0.8 to 2.0. The strength and frequency 
parameters are defined as r4bQ and 2Af/Afj respectively, 



where M is the Mach number of fundamental, b the ncnlin- 
iarity parameter, Q the quality factor of the resonator, 
Af the frequency away from fundamental resonance, and 
A fi the band width at the half-pov/er points . When these 
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results are compared to the theoretical model of Coppens 
and Sanders, it is seen that v;hile the theory accurately 
predicts the magnitude and shape of the harmonic content, 
it consistently underestimates the frequency at which each 
harmonic peaks. In addition, the theory fails to predict, 
except qualitatively, the phase angles. 
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INTRODUCTION 



Coppens and Sanders [1,2] have developed a one- 
dimensional, nonlinear wave equation with dissipative terms 
corresponding to those encountered in a rigid-walled closed 
tube with a large length-to-diameter ratio. An attenuation 
constant of the form 



a = a,//n~ 
n 1 



is assumed, where is the absorption coefficient associa- 
ted with the frequency nw which characterizes the V7all 
losses , R is the radius of the tube and L is the length of 
the tube. With this additional loss mechanism the v;ave 
equation takes the form 
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V = kinematic viscosity, 

Pr = Prandtl number, 

u = angular frequency of the fundamental, 
n = harmonic number. 
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= the thermodynamic speed of sound, 
a-j^ = absorption coefficient for the lowest frequency 

and 

U = particle velocity = + Z U . 

n 



To arrive at the above equation it was necessary to make the 
simplifying assumptions that the wall losses were much 
greater than the bulk losses and that the Mach number 



M = 



P,/P c 
r o 



2 



is much less than unity. is the peak value of the funda- 

mental of the' pressure at the rigid end of the tube. 

A Fourier series solution of the form 



P(X;,t) CO „ 

= — - M S P. (n) Cos In'?! (1- Sin (n(ot+<t>n ) 

p C ^ n=l *- 

o o 

is assumed with the tube excited with a pressure distribu- 
tion given by 

2 

Pi (Xit) = J^PqC^ Cosit (1-x/L) Slnwt. 



When this solution was substituted directly into the non- 
linear wave equation, the following set of algebraic 
equations resulted: 

H^r(Cos0^)S^ - (Sine^)cl 
n L n n n nj 

, n-1 CO 

= Z (S .C.+C .S.)- Z (S„ , -C .S . ) } 



2t '2 . , 
s D = 1 



n-D 3 n-3 3 



j=l 



n+3 j n+: 3 



and 
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where 



and 






r, T n-1 

{i I (C .C. - S .S.) 



- E (C , .C . + S ^ .S .)} , 



j = l 



n+D 3 n +3 3 



n 



~ R Sin <i> r C = R Cos <b , 
n n n 



a = 2/sp, 



6 = Tan 

n 



-1 






1 - n - 2Aa3/o3^6 



n ^tna^/a^ 



H 



n 



a. 



03 = 

r 

SP = 



Tan“^ [A/B] 



, ^ 2 . 

(A"+B ) , 

infinitesimal amplitude absorption coefficient 
in free space, 

angular frequency at resonance , 
strength parameter 



Mbk/2 ^ = Pj^bk/(2p^C^ 



SP 

k = wave number 



Aw = w - w 



These equations were then programmed for computer solution. 
For given values of Aw and SP , the computer calculated 
values of At present, solutions have been 

found for strength parameters up to 0.75. 
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Previous investigations in this area have been carried 
out by Sanders [1,2] and Beech [3]. Both investigators 
examined the harmonic content of the pressure waveform at 
the rigid end of the tube near resonance conditions. 

Sanders employed a wave analyzer to obtain the frequency 
response for the first four overtones. Good agreement with 
theory was found but no phase inform.ation was obtained. 
Beech used photographic and graphical techniques to provide 
numerical analysis of pressure waveforms near resonance. 
Phase information, however, was obtained only for cases at 
resonance. Again good agreement with theory was obtained 
for cases approaching the onset of shock. 

The present investigation examines both the harmonic 
content and phase relationships for strength parameters as 
high as 0.75. Again numerical methods were employed. 
Whereas the graphical techniques of Beech were used to 
generate 64 data equally spaced throughout the pressure 
wave form, the technique in this investigation employed 
analog-to-digi tal conversion by computer to provide better 
than 1000 data to represent the pressure waveform. With 
the pressure waveform digitized to such an extent, it was 
possible to obtain amplitude and phase information for all 
harmonics with amplitudes greater than one percent of the 
fundamental . 
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II. 



EXPERIMENTAL CONSIDERATIONS 



A. APPARATUS 

A block diagram of the experimental system used in this 
investigation is shown in Figure 1. Finite-amplitude 
standing waves were generated in air at ambient conditions 
contained within a tube six-feet long with an inner diameter 
of 2.250 in. and a v/all thickness of 1.125 in. The air was 
excited at one end by a piston which was driven by an M-B 
Electronics Model EA 1500 exciter capable of a maxim.um 
no-load acceleration of 124g. The FA 1500 was driven by two 
M-B electronics Model 2120MB power amplifiers operating in 
parallel. The maximum attainable acceleration with the 
weight of the piston as load was about 50g. The frequency 
and power level were controlled by a General Radio Type 
1161-AGC coherent decade frequency synthesizer. With this 
synthesizer it was possible to select frequency to better 
than 0.001 cycles/sec. (Monitoring the frequency v/ith a 
Hewlett-Packard Model 521C frequency counter showed the 
frequency stability to be better than one part in 10 million.) 
The acceleration of the piston was determined by measuring 
the output an Endeveco Model 2215 accelerometer implanted 
within the piston. When the piston was v/ithdrawn from the 
tube and driven in free air at the approximate experimental 
frequency and acceleration, the harmonic distortion v;as 
measured to be less than 0.3 percent. The other end of the 
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tube was capped with a thick plate which was securely bolted 
to the tube. The sound pressure level at the rigid end of 
the tube was detected with a 1/4'^in. diameter Bruel and 
Kjaer Type 4136 condenser microphone mounted in the cap so 
that its diaphragm was flush with the end of the tube. The 
microphone response was flat to within +0.5 dB from 80 Hz 
to 20 kHz. The outputs of the microphone and accelerometer 
were monitored with Hev7lett-Packard 400D vacuum-tube volt- 
meters. The microphone output also went to a Hewlett-Packard 
Model 302A wave analyzer so that the first tv;o overtones 
could be monitored. The microphone output was recorded by 
a PI 6200 General Purpose Portable Instrumentation Tape 
Recorder. Recordings were made in the FM mode at a tape 
speed of 37.5 in/s so that, the frequency ' response of the 
tape recorder was + 1 dB from DC to 10 kHz. Although the 
signal was recorded at 37.5 in/s it was necessary when 
digitizing the signal to play back at a speed of 3.75 in/s. 

To insure that this change in tape speeds did not alter 
the recorded frequencies, a test signal was recorded at the 
higher speed and then played back at the slow speed. The 
two signals were observed on a dual trace oscilliscope 
where it was observed that the frequency of the test signal 
was 10.0 times that of the signal played back at 3.75 in. 
per second with no discernable drift between the two signals. 

B. MICROPHONE CALIBRATION 

Investigation of the standing waves at various strength 
parameters requires an absolute calibration of the microphone 
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system. The microphone system is defined here as the micro- 
phone, preamplifier and the equipment loaded onto the output 
of the preamplifier. The system sensitivity is defined as 



where is the peak voltage output of the preamplifier and 
P ts the peak pressure at the microphone. Two methods v/ere 
employed to determine S^. The first and simplest used a 
Bruel and Kjaer Model 4220 pistonphone which produces a 
known sound pressure level at 250 Hz. The sensitivity v/as 
found to be = (-77.1 + 0.2)dB re 1 volt/ybar, where the 
uncertainty is that inherent in the pistonphone specifications. 
The second method makes use of theoretical predictions by 
Coppens for small strength parameters as verified by Beech 
[3] . Coppens shows that at resonance the ratio of the peak 
pressure of the second harmonic to the square of the peak 
pressure of the fundamental is given by 



\Vhen the definition of sensitivity is substituted into the 
above equation the result is 



was then found by exciting the tube at resonance and 
measuring the fundamental and. second-harmonic voltages . By 
this method = (-77.2 + 0.3)dB re 1 volt/ybar where the 
uncertainty is due largely to the error in determining the 
attenuation constant. 



S = V /P 
m m' 




=!» = (e<»v2)/(4/2 p^c^cv^). 
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C. MICROPHONE OUTPUT AND STRENGTH PARAMETER 

In this investigation, where the pressure waveform was 
to be analyzed at various strength parameters, it was 
necessary to find a relationship between the strength 
parameter and the observable quantities. Using the 
equations 



S SPL + 21 log 



m 



M = 



V 

rTTs 

p 



2 P 2 

p C o pc 

^o o o o 



+ 72.0, 



= SPL - 197.0, 



and 



M = SP(2a/K) 
b 



one can show that the rms value of the fundamental of the 
microphone output (v) is related to strength parameter and 
the attenuation constant by 



V = 136 . 0 a • SP . 



The constant 136.0 is relatively insensitive to small 
temperature variations about room temperature. 

D. ATTENUATION CONSTANT 

The attenuation constant for infinitesimal amplitudes 
has been shown by Beech to be 




where Af is the frequency difference betv/een the half power 
points on the resonance curve (taken for constant accelera- 
tion amplitude of the piston) . 
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E . FREQUENCY PARAMETER 



The frequency parameter is defined by 



where 



and 



FP 2Aw/w^6j^, 



Aw = CO - w 






This frequency parameter is related to the experimental 
frequency increment (Af ) by 

Af •= Af • FP/2, 

G X 

where Af is the frequency difference between the half power 
points of the infinitesimal resonance curve as used in 
calculating in this investigation Af was found to be 
1.15 + O.OIH^ so that 

Af = 0.58 FP. 
ex 
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III. 



COMPUTER ANALYSIS 



A . SAMPLE RATE 

Many computer programs exist which, when supplied with 
data points defining a periodic waveform, will give the 
Fourier coefficients. The first consideration, generally, 
is the number of data required per cycle. Theoretically 
this number should be at least twice the number of harmonics 
desired. However, the question of accuracy of the computer 
ptedicttons of harmonic amplitudes arises immediately. As 
might be expected, the accuracy of the harmonics determined 
tails off faster as fewer samples per cycle are used. This 
then leads to the problem of first dctcrm.ining hov; many 
significant harmonics are desired. For finite-amplitude 
standing v/aves of large strength parameter the waveform 
qualitatively approaches a function whose harmonics content 
falls off inversely with harmonic number. If such a vjave- 
form is passed through amplifiers and storage devices before 
being analyzed, harmonics below a given amplitude will be 
lost in the inherent noise. The device in this investiga- 
tion which was most seriously restrictive in this respect 
was the PI 6200 FM tape recorder. It was found that the 
dynamic range of this instrument was about 42 db which 
restricted the investigation to harmonics with amplitudes 
greater than about one percent of the fundamental. For a 
periodic ramp function this restriction occurs at about the 
30th harmonic. 
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Determination of the number of samples per cycle to be 
used was based on the arbitrary criterion that the error 
in the 30th harmonic be less than 0.1 percent for an ideal 
ramp function. Applying several computer algorithms for 
determining Fourier coefficients to an ideal ramp function 
one finds that 1000 samples per cycle are required for the 
above accuracy. 

B. THE DIGITAL-COMPUTER PROGRAM 

The primary function of this program was to generate 
normalized Fourier coefficients for a periodic function 
given digitized data on a seven-track tape. Secondary to 
this, the program changed the form of the coefficients to 
A(n) Sin (Nwt + (|) (n) and computed A(n)/A(l) for each harmonic. 
Optional functions provided by the program were the averag- 
ing of several periods of data to reduce the affects of 
random noise. An optional plot of the input data was avail- 
able along with points reconstructed from the calculated 
Fourier coefficients. This plot provided a visual check on 
the input data and the computed coefficients. 

The digitized data provided by the SDS.9300 were written 
on seven-track tape. However, the binary representation of 
numbers on the seven-track tape differs from those on the 
standard nine-track tape used by the IBM 360 system. To 
accomodate the seven-track tape, a conversion program 
(subroutine FORM) was used to convert the data. The con- 
verted data were not floating point but integers and required 
further conversion with the use of a scaling factor. 
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Prime consideration was given to the types of algorithms 



commonly used in generating Fourier coefficients. Two of 
the more widely used methods are the Cooley-Tukey [5] Fast 
Fourier transform and the recursive technique of Goertzel 
[4] . Both methods were tested for accuracy in predicting 
the coefficients of known waveforms. The two methods were 
found to be equal. One advantage of the Cooey-Tukey method 
is its minimal use of computer time. Its computation time 
is about an order of magnitude less than that required by 
the Goertzel method. However, when this time saving is 
compared with the overall required computer time for analysis 
the advantage diminishes. The final consideration and, in 
the end the most important, was the flexibility in using 
data. The Cooley-Tukey method required an integer pov?er-of- 
two number of points. To be able to satisfy this condition 
along with the requirement for 1000 or more samples per cycle 
complicated the digitizing process. The Goertzel method 
required only an odd number of samples and for this reason 
was used in the computer program. 



C. THE ANALOG TO DIGITAL CONVERSION (A-D) • 

The conversion of the analog data was accomplished with 
the use of a hybrid computer consisting of the SDS 9300 
high-speed digital computer and the Ci 5000 general-purpose 
analog computer. Use of the system for A-D conversion 
requires a FORTRAN program for the SDS 9300 and properly 
constructed analog and logical patch boards for the Ci 5000. 
In addition, several peripheral equipments are required. 
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These are the card reader, magnetic tape transport, line 
printer and, if desired, a pen recorder. 

Aside from the question of sampling, which was deter- 
mined earlier, the important question in the A-D conversion 
is accuracy. The clock in the computer is accurate to one 
part in 100 million, well within any accuracy requirements. 
The A-D converter is designed to accept signals between 
+ 100 volts and divide this interval into 16,384 parts 
providing a resolution of 12.2 mV. For a properly amplified 
analog signal, the dynamic range of the A-D converter is 
much greater than the PI 6200 tape recorder. 

The analog patch board (Figure 3) , as used for A-d con- 
version, served to amplify the analog signal prior to con- 
version and distributed the signal to the converter and 
oscilliscope . Additionally, a patch was provided linking 
the Digital-to-Analog converter with a pen recorder. This 
additional function is explained in detail in the pro- 
cedures section. 

The logic patch board (Figure 4), had three functions. 
The first provided a sampling frequency (10,000 Hz) to the 
A-D converter. This frequency represents the approximate 
maximum sampling rate allowed. Since the approximate 
X'esonAnt frequency of the tube was 94 Hz it was necessary to 
play back signals from the PI 6200 at one tenth the recorded 
speeds to acquire the desired number of samples per cycle. 
The second function provided a triggering voltage through 
an external switch (DS-1) which was used to start the data 
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conversion process. Associated with this function is an 
externally controlled delay (DFOO) which insured that the 
triggering voltage remained on long enough for one set of 
data (one record) to be written on the seven-track tape. 

The last function provided another triggering voltage 
through an external switch (DS-2) . When the last record was 
to be written on tape switch DS-2 was turned on while 
switch DS-2 v/as triggered. 

D . SYSTEM TEST 

Prior to analyzing experimental data, it was necessary 
to test the entire system process of A-D conversion and 
computer analysis. The results of this test for a ramp 
function is shown in Figure 5. 
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IV. 



PROCFDURFS 



A. PRERUN PROCEDURES 

Beech pointed out that the resonant frequency of the 
tube changed by 0.16 Hz for each degree centigrade of 
temperature change. To minimize this affect, all data 
gathering runs were made during the evening hours v/hen 
the room temperature remained relatively constant for 
several hours. Typical resonant frequency variations during 
these hours were of the order of 0.01 Hz per hour. Several 
hours prior to taking data the equix^ment was turned on 
and allowed to stabilize. During this warmup period, the 
PI 6200 tape recorder was calibrated to insure that the 
FM system was properly aligned and to insure that the 
signals being recorded were not overdriving the tape recor- 
der preamplifier and thereby introducing unwanted signal 
distortion. Next, the alignment of the piston in the free 
end of the tube was checked by observing the output of the 
accelerometer on an oscilloscope and measuring, with a 
wave analyzer, the second- and third-harmonic content. If 
necessary, the harmonic content could be reduced by carefully 
adjusting the alignment of the piston so that the distortion 
from any harmonic was less than 2 percent. To minimize 
leakage of sound from the tube, the piston was fitted with 
an "O" ring and lubricated. The harmonic distortion intro- 
duced by the piston was sensitive to the type of lubricant 
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used. It was found that a light turbine oil minimized the 
distortion while providing an adequate seal. 

Just prior to taking each set of data, it was necessary 
to determine the resonant frequency and the attenuation 
constant. For this measurement the piston was excited at 
a constant value of acceleration low enough so that the 
strength parameter at resonance was less than 0.03. From 
the frequency difference between the half power points, the 
attenuation constant was determined. The resonance frequency 
was taken to lie half way between the half power points. 

B. DATA TAKING PROCEDURE 

Corresponding to the desired strength parameter was 
a given value of the pressure fundamental as shown earlier. 
With the piston excitation set to the proper level to keep 
the pressure fundamental constant, data v;ere recorded at 
0.1 Hz intervals for twenty frequencies from 0.6 Hz below 
resonance to 1.3 Hz above. At each frequency, the signals 
from the microphone and accelerometer were recorded on the 
PI 6200 at high speed (37.5 in/s) for about 10 seconds 
(30 ft.), and the second-harmonic amplitude of the pressure 
was noted to be later compared with the computer analysis. 
Since the taking of an entire set of data at one strength 
parameter required from 30 minutes to an hour, at the end 
of the run it was necessary to check the resonant frequency 
and attenuation constant to insure that the variation in 
resonant frequency was substantially less than one frequency 
interval (0.1 Hz) and that the change in attenuation constant 
was less tlian one percent. 
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The theoretical predictions (8) are compared with the 
experimental results in Figures 6 through 13 for each value 
of strength parameter, plots are shown of the harmonic 
amplitude vs. frequency parameter and of the phase difference 
between consecutive harmonics vs. inverse harmonic number. 
Figure 6 shows the actual data obtained for three different 
experimental runs. The degree of scatter is a demonstra- 
tion of the excellent reproducibility of the results. Both 
experimental and theoretical results are limited to harmonics 
with amplitude greater than 1% of the fundamental. 

In comparing the experimental values of the harmonic 
amplitudes with the theoretical predictions for a strength 
parameter of 0.25 it is seen that the agreement is fairly 
good but that there does seem to be a slight systemat c 
difference. The corresponding curves for higher strength 
parameters display this difference to a more marked degree. 

It is seen that the experimental results peak at about the 
predicted values of that the shape of the 

curve for each harmonic is about as predicted. However, 
the experimental results peak at a value of frequency 
parameter significantly greater than predicted and this 
difference increases with harmonic number. It is interest- 
ing to note that at the frequency of maximum distortion of 
the pressure waveform (at approximately FP = 0.5) theory 
accurately predicts the harmonic amplitudes . 

The comparison between the experimental and theoretical 
phase angles displays an even greater discrepancy. V^hile 
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values of do seem to be linearly related to 1/n , 

the values for the second harmonic are greatly different 
and the slopes of the experimental curves are significantly 
greater than predicted. 

In conclusion, it appears that while the present theo- 
retical model is able to predict the magnitude and shape of 
the harmonic content of the distorted waveform, it consis- 
tently underestimates the frequency at which each harmonic 
peaks. This underestimatioji increases with harmonic number 
and strength parameter. In addition, it fails to predict, 
except qualitatively, the phase relations in the distorted 
wave . 



22 



APPENDIX A 



THE A-D PROGRAM PARAMETERS AND OPTIONS 

The purpose of this program was to convert up to four 
channels of analog data into a series of digital points on 
seven-track magnetic tape. The program had several input 
parameters which had to be provided to the computer via 
the teleprinter after the FORTRAN program had been read into 
the system. 

NSAMP (input 1.) was an integer which corresponded to the 
number of data or samples to be taken from each channel. 

For this investigation a sampling rate of 10 kHz was used 
and NSAJ1P v;as chosen to span one cycle of the input waveform. 

NCHAN (input 2.) was an integer which specified the 
number of channels of analog data being converted. 

NREC (input 2.) specified how many records were to be 
written on tape each time sampling was initiated. A record 
was defined as the data set corresponding to the product of 
NCHAN and NSAMP. Because of the high sampling rate used it 
was necessary to set this value to 1. 

ITAPE (input 4.) designated which magnetic-tape trans- 
port held the seven-track tape. The integer used here had 
to correspond that selected on the "unit select" switch on 
the particular tape unit used. 

NDEL (input 5.) determined how rapidly the data on the 
seven-track tape was read back when using program option 7 . 
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An example of the input parameters as typed on the 

I 

teleprinter might be 

NREC = 1, NSAMP = 1000, NCHAN = 3, ITAPE = 1, 

NDEL = 20* (carriage return) . 

The control symbol * indicated to the computer that the 
input parameters had been set. Care was taken when choos- 
ing NSAMP and NCHAN to insure that their product did not 
exceed the first dimension of IBUF found in the first FORTRAN 
program statement. 

When the parameter list had been given the teleprinter 
responded with 

"OPTION = (II) " . 

The computer v/as then given a one-digit number corresponding 
to the list of available options. 

RESTART (option 1.) was used whenever any one or all of 
the input parameters was to be changed. 

GO (option 2.) was used when analog to digital conver- 
sion was to be done. Upon receiving this option the compu- 
ter awaited the initiating trigger from switch DS-1 to 
start the actual conversion process. Each trigger from DS-1 
caused one record of data to be written on tape. When a 
sufficient number of records were written, DS-1 was 
triggered with DS-2 in the "up" position. This caused the 
last record to be written followed by a teleprinter response 

"OPTION = (II)". 
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END OF FILE (option 3.) was used to put a mark on the 
tape to separate or isolate groups of records. Groups of 
records separated by END OF FILE marks are called files. 

SKIPFILES (option 5.) were used when files already on 
tape were to be skipped over. When used, a four-digit 
number was given the computer to specify the number of 
files to be skipped. 



LIST ONE RECORD (option 6 . ) was used in reading data 
from the seven-track tape. Use of this option caused data 
written on tape to be written on the page printer after 
having specified the number of data to print. 

PLOT ONE RECORD (option 7.), when called, caused one 
record to be read from the tape, converted to an analog 
signal, and fed to terminals T420 through T423 on the analog 
board (corresponding to channels 1 through 4). The signal 
could then be monitored on various recording devices. 

Because of the short duration of the signal it was best 
displayed on the available pen recorder whose input terminals 
are found on the analog board (Figure 3) . The time base 
of the signal returned from the tape could be lengthened or 
shortened by varying the value of the input parameter NDFL . 
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Figure 1. Apparatus. 
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Figure 2. PI 6200 Dynamic Range. 



27 




28 



Figure 3. Analog Patch Board. 
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Figure 4. Logic Patch Board. 
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Figure 6. Point Scatter 
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Figure 7. Harmonic Amplitudes SP = 0.25. 
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Figure 8. Harmonic Amplitudes SP = 0.50. 
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Figure 9. Phase Plot SP = 0 . 50 . 
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Figure 10. Harmonic Amplitudes SP = 0.75. 
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Figurell . Phase Plot SP = 0.75. 
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Figure 12. Harmonic Amplitudes SP 
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Figure 13. Phase Plot SP = 1.00 
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